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Path integral techniques for the pricing of financial options are mostly based on models that can 
be recast in terms of a Fokker-Planck differential equation and that, consequently, neglect jumps and 
only describe drift and diffusion. We present a method to adapt formulas for both the path-integral 
propagators and the option prices themselves, so that jump processes are taken into account in 
conjunction with the usual drift and diffusion terms. In particular, we focus on stochastic volatility 
models, such as the exponential Vasicek model, and extend the pricing formulas and propagator 
of this model to incorporate jump diffusion with a given jump size distribution. This model is of 
importance to include non-Gaussian fluctuations beyond the Black-Scholes model, and moreover 
yields a lognormal distribution of the volatilities, in agreement with results from superstatistical 
analysis. The results obtained in the present formalism are checked with Monte Carlo simulations. 



I. INTRODUCTION 

It is well known that the pioneering option pricing 
theory of Black and Scholes [l| and Merton Q fails 
to reflect some important empirical phenomena. Many 
studies have been conducted to modify and improve the 
Black-Scholes model. Among others, popular models in- 
clude, (a) the local volatility models Q; (b) the stochastic 
volatility (SV) models [j-H ; (c) the SV and stochastic in- 
terest rate models |7H10l| ; (d) the jump diffusion models 
pdT[l3j ; (e) models based on Lev y p rocess [T3 - fl8| : and 
(f) the SV jump diffusion models HM! • 

Inspired by [15|, Il9l42l[ |26| we will focus on the latter 
class of models. For example, Cont and Tankov [l5[ and 
Gatheral [26| motivate that the combination of jumps in 
returns and SV makes it possible to calibrate the implied 
volatility surface, without using time dependent param- 
eters. Jumps make it possible to reproduce strong skews 
and smiles at short maturities while SV provides for the 
calibration of the term structure, especially for long-term 
smiles. 

In this article we will present a method that makes it 
possible to extend the Fourier space propagator of a gen- 
eral SV model to the Fourier space propagator of that SV 
model where an arbitrary jump process has been added 
to the asset price dynamics. Thereby we contribute to 
the existing work on Fourier transform methods applied 
to option pricing. For example in fTTI j jump diffusions are 
treated and prices for some exotic options are obtained. 
In 27] the Heston model is extended with a jump process 
for the asset price. In [25[ the Heston model is extended 
with arbitrary jump processes in both the asset price and 
the volatility process. 

As an application, we investigate a model where we 
assume that the stochastic volatility follows an exponen- 
tial Vasicek model [28|,[29j]. To the best of our knowledge, 
for this model no closed form formulas for the propaga- 



tor or the vanilla option price exist yet. Making use of 
path integral methods [9l. KM l3l| we derive approximative 
closed form formulas for the propagator and for vanilla 
option prices for this model (for more information about 
methods from physics applied to finance see for exam- 
ple [32T - f34j ). Using Monte Carlo (MC) simulations we 
specify parameter ranges for which the approximation is 
valid. Using the above mentioned method we extend the 
propagator of this model to the propagator of this model 
extended with jumps in the asset price which leads also 
to closed form pricing formulas in this extended model. 
Also these last results are checked with MC simulations. 

This paper is organized as follows. In section |TT] we 
present the method for extending the propagator of a 
general SV model to the propagator of that model with 
jumps in the asset price. In section lLTIl we present an ap- 
proximative propagator for jump diffusion models where 
the volatility is assumed to follow an exponential Vasicek 
model. Section HVl is devoted to European vanilla option 
pricing, as well as comparisons with MC simulations. In 
this section we also give parameter ranges for the ap- 
proximation made in the exponential Vasicek model to 
be valid. And finally a conclusion is given in section fVl 



II. GENERAL PROPAGATOR FORMULAS 

A. Arbitrary SV models 

We assume that the asset price process S(t) follows the 
Black-Scholes stochastic differential equation (SDE): 



dS(t) = rS(t)dt + (j{t)S{i)dWx{t), 



(1) 



in which r is the constant interest rate and the volatility 
a(t) is behaving stochastically over time, following an 



2 



arbitrary stochastic process: 

da{t) = A(t, o(t))dt + B(t, a(t))dW 2 (t). 



(2) 



Here and in the rest of the article Wj = {Wj(t),t > 
0} (j — 1,2) are two correlated Wiener processes such 
that Cov[dWi(t) dW 2 (t) ] = pdt. 

Eg .([!]) is commonly expressed as a function of the lo- 
greturn x(t) = \nS(t), which leads to a new SDE: 



dx(t) 



a 2 (t) dt + a(t)dWi(t). 



(3) 



To deal with the pricing problem, we need to solve for 
the propagator of the joint dynamics of x(t) and <j(t). 
The propagator, denoted by V(xt, ctt, T\ xq, ao, 0), de- 
scribes the probability that x has the value xt and a has 
the value ot at a later time T given the initial values xq 
and do respectively at time 0. It satisfies the following 
Kolmogoroff forward equation: 



dV 
df 



8 
dxr 
d 



<r--a 2 T )V 



1 d 2 



2dx\ 



d 2 



^ 8xt d(JT 

with initial condition 



[a T B(T, a T )V\ 



(4) 



V(x T ,Cr T , 0| X , CTq, 0) = S(x T - X ) 5((J T - <T ). (5) 



B. SV jump diffusion models 

A general SV jump diffusion model is obtained by 
adding an arbitrary jump process into the asset price 
process (see for instance |2Q|). That is, equation (fT]) be- 
comes 

dS{t) = fiS(t)dt + a{t)S(t)dWi(jt) + (e J - l) S{t)dN(t), 

(6) 

where N — {N(t),t > 0} is an independent Poisson pro- 
cess with intensity parameter A > 0, i.e. E,[N(t)} = Xt. 
The random variable J with probability density ■uj( J) 
describes the magnitude of the jump when it occurs. 

Here the risk- neutral drift \i = r — A m? is no longer 
the constant interest rate r, rather it is adjusted by a 
compensator term Am- 7 , with m J the expectation value 
ofe J -l: 



m? = E \e J - ll = 



+ 00 



(e J - l)m(J)dJ, (7) 



so that the asset price process constitutes a martingale 
under the risk neutral measure. And the logreturn x(t) 
follows a new SDE: 

dx{t) = (r- Xm j - l -a 2 {t)) dt + a(t)dWi{t) + JdN(t). 



(8) 



Given the same arbitrary SV process ([2]), the new prop- 
agator of this model, denoted by Vj(xt, <Jt, T\xq, <to, 0), 
satisfies the new Kolmogoroff forward equation (see for 
instance (35[) 



9Vj 
dT 



d 
dxr 
1 d 2 



2dx 2 T 
1 d 2 



Bar 



-A(T,a T )Pj 



2 da 2 
d 2 



[B 2 (T 7 a T )V.i] 



+P- 



[a T B{T,a T )Vj 



dxr 9ot 

/+00 
[Vj(x T - J) - Vj(x T )} w(J)dJ.(9) 
-00 

If we write the propagator of the arbitrary SV model 
as a Fourier integral (here and below, i is the imaginary 
unit) 

V(x T ,a T ,T\x ,ao,0) 
+°° dp 



2vr 



W*T-*o) F fa,a ,r, Pi T), (10) 



then the propagator of arbitrary SV jump diffusion mod- 
els can be written as 

Vj(x T ,a T ,T\x o ,ao,0) 
30 dp 



2tt 



e ip ^ T - x ^F(a T ,a ,r,p,T)e u ^' T \(ll) 



where 



U(p,T) = \T 



-\-oo 



-ipj 



+ ip (e J - 1)] w{J)dJ. 

The proof of this statement is given in the Appendix 
„ Note the relation between propagators (fTU|) and pTT) . 
The only difference between them is the factor e u ^ p ' T \ 

If this is applied to the propagator of the Heston model 
[9|, the propagator of the Heston model with jumps is 
obtained. This propagator is similar as the one derived in 
Ref. [25(. Furthermore the above described method can 
be combined with the method described in Ref. 9] for 
finding the propagator of a model including both SV and 
stochastic interest rate. In particular extending the result 
of Ref. [9J for the Heston model with stochastic interest 
rate to include jumps again only involves multiplying the 
propagator with e u ^ p ' T ^ as in pT|) . In the next section, as 
an example of the method of this section the volatility of 
the asset price will be assumed to follow an exponential 
Vasicek model. 



III. EXPONENTIAL VASICEK SV MODEL 
WITH PRICE JUMPS 

The Heston model assumes that the squared volatility 
follows a CIR process which has a gamma distribution 
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as stationary distribution. This assumption should be 
compared with market data. Attempts to reconstruct the 
stationary probability distribution of volatility from the 
time series data (among others, see Refs. [H, [29|, 1331 ]) 
generally agree that the central part of the stationary 
volatility distribution is better described by a lognormal 
distribution. 

Due to the different structure in path-behavior be- 
tween different models, Schoutens, Simons and Tistaert 
find that the resulting exotic prices can vary significantly 
17]. So an investigation into an alternative model which 
fits market data better is meaningful. 

Furthermore the model will serve here both to demon- 
strate the use of path integral methods in finance and to 
illustrate the method of section [Til 

When a(t) is assumed to be an exponential Vasicek 
process (used for example by Chesney and Scott [36|). 
this results in the following two SDEs 



where the Lagrangians are given by: 



dS = rSdt + aSdWi, 



(13) 



1 



da = cr ( [a - In a] + ^p 1 ) dt + 1 adW 2 . (14) 



This model has a lognormal stationary volatility distri- 
bution and we will denote it by the LN model, the prop- 
agator for this model will be denoted by Vln- In this 
model In a(t) is a mean reverting process, with the 
spring constant of the force that attracts the logarithm 
of asset volatility to its mean reversion level a. Again 7 
is the volatility of the asset volatility. As far as we know, 
there is no closed form option pricing formula for this 
model. In this section, we will give an approximation 
for the propagator of this model. In the next section we 
will give an approximation for the vanilla option price 
and determine a parameter range for which the approxi- 
mation is good. The derivation starts with the following 
substitutions: 



y(t) = X (t) - - e 2(t) - rt, 
1 

z(t) = In a(t), 



(15) 
(16) 



where x(t) is defined as before. This leads to two uncor- 
rclatcd equations: 



dy 



dz 



1 



n 7 ( 0(a — z) 7" 



dt 



+eVl-M, 
0(a- z)dt + jdB 2 , 



(17) 
(18) 



where B\ and B 2 are two uncorrelated Wiener processes. 
Since these equations are uncorrelated, the propagator 
'PlnIut, zt \ yo, zq) is given by the following path integral 

V L N(yT,z T \ yo,zo) 
= I P: [ / Vye-fo^M dt J e -/o T *H* (19) 



C[z] 



y + 



2(1 -p 2 )e 2z 

{z~0(d-z)} 2 
2^2 2' 



(20) 
(21) 



The first step in the evaluation of (| is the integration 
over all y paths. Because the action is quadratic in y, 
this path integration can be done analytically and yields 

VLN(yT,z T \yQ,zo) 

Vz e~ S « c[z]dt 



x e 



yj2n(l-p*)f*e**dt 
^+p£(f^+|>H! 



(22) 



Note that the probability to arrive in (yT,zx) only de- 
pends on the average value of the volatility along the 
path z(t), in agreement with Ref. [36T ]. With the help of 
a Fourier transform, we rewrite the preceding expression 
as follows 



"P LN (y T ,z T \y ,z ) 
+ °° dp 



2tt 



(l-P )P ~ip fT 2 



x e 



e^dt+i PP j t ( ^ +i y dt (23) 



If ((t) — z(t) — a, then ((t) is close to zero because z(t) 
is a mean reverting process with mean reversion level a, 
This motivates the approximation » 1 + ( + This 
type of approximation is akin to expanding the path in- 
tegral around the saddle point up to second order in the 
fluctuations, as in the Nozieres-Schmitt-Rink formalism 
extended to path-integration by Sa de Melo, Ran- 
deria and Engelbrecht [38]. Now we can work out the 
remaining path integral in (|23[) 



vc - Io T { ^ ~§ + ±e 2 < + B K e<-^e<}dt 

[( CT + ^M) 2 _( Co + ^M) 2 ]_,( 4 _ c g) 



,2 J„ll 



xe 





2 







7T7 2 (1 - e~ 2wT ) 



^ iT+ £My^ Q+1 2My^ 2 
(24) 



-■2, 1 .,.--2uT 



4 



where 



.4 
B 

M 




(25) 
(26) 

(27) 

(28) 



We see that also the integral over the final value Cr can 
be done, yielding the marginal probability distribution: 



Vln(xt\x , Co) 

+°° dp 



ip[XT—XQ—r 



2tt 



T]+B(e c °-l) 



xe 



,7 2 [2u+(3-^ + B 7 2)(l- 



1 + ^ 



UJ 



£ 7 2 



(29) 



where 



2B 1 2 N + u(N - 



7 2 M 



N 



+ (1 



7 2 M 



-2uiT\ 



UJ- 

gy _ bj 4 m 

2 w 

08 + B 7 2 ) 7 4 M 



) 2 - 08 + ^2)^ 



2w 3 



(Co 



7 2 M 



(30) 



(31) 



The goodness of this approximative propagator needs 
the support from MC simulations because of the lack of 
a closed form solution. 

Figure Q] shows the propagators as a function of xt — xq, 
i.e., ln-^ 1 . The full curves come from expression (|2U|) . 
while the marked ones are MC simulation results, with 
time to maturity ranging from three months to five years, 
and correlation coefficients 0, —0.5 and 0.5 respectively. 
Here and in the rest of the article we will set oo equal to 
the long time average of the volatility: 

<j = lim E[<r(t)] =exp(a+^l, (32) 
t^oo [ 4pJ 

which seems a reasonable choice. For these MC simula- 
tions 5,000,000 sample paths are used. 

It is seen that our analytical results fit the MC simu- 
lations quite well. Actually, using the parameters of Fig. 
[TJ and putting expression (|2"9")l for those three cases into 
the left hand side of the Kolmogorov backward equation: 



_ &P 
~ df 
dV 



r -l p 2(Co+a) 
2 



dV 

dx 



2 dxl 



1 2 d 2 V 

-7 

2 ' 



pe Co+a 7 



d 2 v 

9x Co 



= 0, 



(33) 




FIG. 1: Propagator V(xt\xo, Co) as a function of xt — xq. The 
full curves are our analytical results, while the symbols repre- 
sent Monte Carlo simulations. T — 0.25y, p = (crosses). T 
= ly, p = -0.5 (circles). T = 5y, p = 0.5 (triangles). For the 
other parameters the following values are used for the three 
figures: /3 = 5, 5 = —1.6, 7 = 0.5, r = 0.015. 



we find that, for different Xt values, the absolute values 
are all in the order of 10 ~ 7 or even smaller. In section 
HVBI we come back to the discussion concerning the good- 
ness of our approximation. 

According to the discussion of Section [Til an extension 
of this model to the one with price jumps is straightfor- 
ward: the new marginal probability distribution would 
be: 



V LNJ {x T \0Xo) 



-Hoc 



dp ip\x T ~x Q ~rT]+B(e c Q - l) 



00 



2tt 



2 \ 2 



gC 2_^ C() + ^M) 



xe 2 ~i 

e ^2[ 2 „ + ( 3 _ Iil + B 7 2 )(1 _ c -2 IJ T ) ] 



+ \ H-^-A + B^ +J ^_l T 



X- 



1 



1-e 



-2uT 



\P-u + B~ f 2 



xe 



XT /-~[ e ~ ,pJ -l+ip(e J -l)]w(J)dJ 



(34) 



where the same notations as in Eq.(f2U)) are used. 
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IV. EUROPEAN VANILLA OPTION PRICING corresponding U{p,T)'s: 



A. General Pricing Formulas 

If we denote the general marginal propagator by 

" + OO 



,U(p,T) 



T(x T \x ,a )= J °° ^.^-o-rT) F ^ T)[ 

^ . (35) 

then the option pricing formula of a vanilla call option C 
with expiration date T and strike price K is given by the 
discounted expectation value of the payoff: 



C = 



(e XT - K),V(x T \x ,a )dx T 



+ i 



2tt 



(36) 



where 



G(p) = S F(p + i,T)e u ^ +i ^ 
-Ke- rT F( P) T)e u ^ T \ 



(37) 



and (x) + — max(x,0). 

Here we have followed the derivation outlined in Ref. 
[39| . In particular for the LN model F(p, T) equals: 



F{p,T) = 



2-^ 



S ( e <0 l)+ 72[2w+(/3 _^ + s 7 2) (1 _ c -2u,T ) ] 

x -.(38) 



1 



l-e- 



2^ 



\P-U + B 7 2 ] 



At this stage one needs to specify the PDF for the 
jump sizes. Merton [l3| and Kou [HI proposed a normal 
distributed jump size, denoted by wm{J), and a asym- 
metric double exponential distributed one, denoted by 
vjk(J), respectively: 



TDM (J) 
ZU K (J) 



1 



p + —e~^ J 0(J) 
V+ 

1 — j 

+ p- — e*>- 6(-J) 
77_ 



(39) 



(40) 



For the Merton model v is the mean jump size and S is 
the standard deviation of the jump size. For Kou's model 
< 77+ < 1, rj- > are means of positive and negative 
jumps respectively. p+ and p- represent the probabilities 
of positive and negative jumps, p + > 0, p_ > 0, p + + 
P- = 1 and 9 is the Heaviside function. 

According to expression (JTSJ) , it is easy to derive their 



Um(p,T) 



XT 



-ipv—^S p 



ip (e"+i* 2 - l) 



U k (jp,T) = XT 



-ip 
P+ 



P- 



1 



- ipr\- 

V- 

1 + 7?- 



1 



(41) 



(42) 



Using expression (13"81 and results (|4T|) . (|4^|) in formu- 
las (|36|) . (|37)l allows to find the price of the vanilla call 
option for the exponential Vasicek stochastic volatility 
with price jumps model. 



B. Monte Carlo simulations 

To test our analytical pricing formula for the LN 
model, we focus on the parameters that most strongly 
influence the approximation. To satisfy the assumption 
that quadratic fluctuations around the mean reversion 
level a captures the behavior of the volatility well, the 
mean reversion speed /? and the volatility 7 of asset 
volatility are crucial. 

The substitution r = j 2 t transforms expression 
into 



dz{r) 



- ? 2 



:{T)}dT + dB 2 (T), 



(43) 



showing that it is actually the parameter c = -4r which 
determines whether the approximation will be good. For 
bigger c values the approximation z (t) ~ a will be better. 

As the correlation parameter p controls the skewness 
of spot returns, we will also consider the typical negative 
and positive skewed cases by taking values —0.5, and 
0.5 for this parameter. On the other hand, the constant 
interest rate r and the mean reversion level a do not 
influence the accuracy of the result a lot, and we just 
assume them to be constant values: r — 0.015 and d = 
— 1.6 « In 0.2. These two parameters seem to be quite 
reasonable for the present European options. 

To get an idea of what is a reasonable range for c, 
and since calibration values for the LN model are not 
available, we took calibration values from the literature 
4, 40] for the Heston model and fitted our model to the 
volatility distribution of the Heston model with those 
parameters. For Q we obtained c w 7 and for [4(| c ~ 18. 
Therefore in Table I we used values for /3 and 7 such 
that c ranges from 4.08 up to 25. We calculated prices 
for S* = 100 and K = 90, 100 and 110. 

The comparison of our analytical solution with the MC 
solution for a European call option in the LN model as 
shown in Table U suggests that for the above mentioned 
parameter values the relative errors are less than 3% and 
most of the time even less than 1%, which is acceptable 
when we take the typical bid-ask spread for European 
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FIG. 2: The upper figures show European call option prices in the LN model (left), the LN model with Merton's jump (middle) 
and the LN model with Kou's jump (right). The red curves are our analytical results and the black crosses are the Monte Carlo 
simulations. The corresponding lower figures give the relative deviations of our analytical results from the MC simulations in 
the unit of percent. Parameter values So = 100, r = 0.015, T = 1, j3 = 5, a = —1.6, 7 = 0.5, p = —0.5, A = 10, v — —0.01, 
S = 0.03, p+ = 0.3, p_ = 0.7, r)+ = 0.02, r\- = 0.04 are used here. 



options into account. Here each MC simulation runs 
20,000,000 times. 

For the basic LN model we can conclude that we found 
an approximation valid up to 3% for parameter values 
c > 7 (We only checked values of c < 25, but for bigger c 
the approximation will only become better), —0.5 < p < 
0.5, T < 1 and 0.9 < K/S < 1.1. 

Finally we consider the vanilla call option pricing in 
LN model combined with Merton's and Kou's jumps, re- 
spectively. Since the jump process is independent from 
the approximation we made, we do not investigate the 
goodness of our approximation as thoroughly as in the 
basic LN model (assuming that, if it is good there it will 
be good here). Figure [5] illustrates our analytical results 
(curves) and the MC simulations (crosses) , as well as the 
relative errors in the unit of percent. Each MC simula- 
tion runs 300,000,000 times. These results suggest that 
the approximation error is typically less than 2%. And 
due to the fact that whenever the degree of moneyness 
(the ratio of the strike price K to the initial asset price 
So) is relatively high, the average bid-ask spread tends 
to be relatively high for call options [4l|, our analytical 
results can serve as an easy way to get a quick estimate 
that is normally accurate enough for many practical ap- 
plications. 



plied to the Heston model, leads to similar results as 
those obtained in Ref. (25[, which gives us confidence 
in the present treatment. The stationary volatility dis- 
tribution of the Heston model, however, does not corre- 
spond to the observed lognormal distribution [28, ^2t| [33[ 
in the market. The exponential Vasicek model does have 
the lognormal distribution as its stationary distribution. 
Therefore we used this model for the volatility to illus- 
trate the method presented in section [TTJ For this model 
no closed form pricing formulas for the propagator or 
vanilla option prices exist. We first derive approximative 
formulas for the propagator and vanilla option prices for 
this model without jumps, using path integral methods. 
This result was checked with a Monte Carlo simulation, 
proving a parameter range for which the approximation 
is valid. We specified a parameter range for which our 
pricing formulas are accurate to within 3%. They become 
more accurate in the limit \ » 1 where /3 is the mean 
reversion rate and 7 is the volatility of the volatility. Fi- 
nally we extended this result to the case where the asset 
price evolution contains jumps. 



Appendix: Derivation of equations (|12|) . 



V. CONCLUSION 

We presented a method which makes it possible to ex- 
tend the propagator for a general SV model to the propa- 
gator of that SV model extended with an arbitrary jump 
process in the asset price evolution. This procedure, ap- 



The proof starts by assuming that a solution for 
Vj(xt,(Jt,T\xo,<jo,0) of the form ((TTJ exists. Below 
we show that this assumption indeed leads to a so- 
lution, which in turn justifies the assumption. Since 

/- ~ Ige^-xo) 9F(a T ,ao,r,p,T) ^ ^ hand 

side of Eq. Q and the derivative operators and g^— 



have no effect on e u ^ p ' T \ it follows that: 



the right hand side of Eq. (|9]) is expressed as 



d 

1 d 2 



2&4 
1 d 2 



d<7T 



-A(T, a T )Vj] 



[B 2 (T,a T )Vj] 



2 da 2 
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^ dxr d(JT 



[a T B{T,a T )Vj 



ip(xT-Xo) ' 

2tt 



dT 



,u( P . 



^A.l) 



Adding the term XmP -j^r-Vj, which is given by 



/ +OC ^Pe lp{xT - Xo) F(a T , c , r,p, T) , 



U(p,T) 



+oo 



x / (e J - l)w(J)dJ, 



as well as the term A [Pj{xt — J) — 'Pj(xt)] zj(J)dJ, 
which is given by 



d P r ir>(x T -x a ) 9F{(7 T , go, r,p, T) 



,£/(p,T) 



— OO 



2tt 



<9T 



* e ^-o)^ ((TTjCTOjr , jPiT)eC /( P ,T) 
Z7T 



xA / [e- !pJ - 1 + ip(e J - 1)] n7(J)dJ. (A.4) 



This, of course should equal the left hand side of Eq.©, 
which is given by 



d P r iv(x T -x a ) dF(a T , op, r,p, T) u{v T) 
2tt dT 



JijpixT-Xo)^^ J) - ( A.5) 



2tt 



de W) 
~df~ 



(A. 2) Expression (|A.4[) equals (|A.5|) when 



dU(p,T) 
dT 



A 



" „ip(xT— xo 



/-l-oo 
[e- lpJ - 1 + ip (e J - 1)] m(J)dJ, 
-OO 

(A.6) 



/e^-^F(a T ,a ,r,p,r) e 
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from which the result (fT2"|) for t/(p, T) follows. 



x / (e~ tpJ - 1) w{ J)dJ, 
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TABLE I: Comparison of our approximative analytic pricing result and the MC simulation value for the LN model. 



Parameter values 



Relative error 



K 


P 


7 


P 


MC value(a) 


Approx.(b) 


(b- a) /a (%) 


90 


-0.5 


1.2 


7 


15.3947 


15.2533 


-0.9185 








8 


15.2979 


15.1731 


-0.8166 








10 


15.1630 


15.0588 


-0.6869 






0.8 


5 


15.1079 


14.9995 


-0.7180 








6 


15.0307 


14.9337 


-0.6475 








7 


14.9776 


14.8855 


-0.6151 







0.7 


2 


15.2486 


15.1982 


-0.3307 








3 


15.0259 


15.0024 


-0.1564 








4 


14.9190 


14.8992 


-0.1328 






0.5 


1 


15.2061 


15.1576 


-0.3187 








2 


14.9030 


14.8882 


-0.0996 








3 


14.7951 


14.7865 


-0.0577 




0.5 


0.3 


1 


14.6051 


14.5035 


-0.6953 








1.5 


14.5524 


14.4815 


-0.4872 








2 


14.5354 


14.4775 


-0.3986 






0.2 


0.5 


14.6015 


14.5398 


-0.4192 








0.75 


14.5609 


14.5098 


-0.3519 








1 


14.5332 


14.4981 


-0.2418 


100 


-0.5 


1.2 


7 


9.4541 


9.2862 


-1.7762 








8 


9.3720 


9.2197 


-1.6253 








10 


9.2599 


9.1274 


-1.4310 






0.8 


5 


9.1537 


9.0346 


-1.3006 








6 


9.0950 


8.9869 


-1.1886 








7 


9.0557 


8.9534 


-1.1291 







0.7 


2 


9.5394 


9.4975 


-0.4395 








3 


9.2906 


9.2704 


-0.2174 








4 


9.1682 


9.1513 


-0.1840 






0.5 


1 


9.4915 


9.4493 


-0.4480 








2 


9.1445 


9.1328 


-0.1285 








3 


9.0235 


9.0155 


-0.0887 




0.5 


0.3 


1 


9.0168 


8.9081 


-1.2051 








1.5 


8.9302 


8.8522 


-0.8737 








2 


8.8886 


8.8246 


-0.7195 






0.2 


0.5 


8.9655 


8.9023 


-0.7048 








0.75 


8.9039 


8.8509 


-0.5955 








1 


8.8628 


8.8247 


-0.4295 


110 


-0.5 


1.2 


7 


5.2749 


5.1365 


-2.6237 








8 


5.2209 


5.0916 


-2.4758 








10 


5.1507 


5.0335 


-2.2756 






0.8 


5 


5.0170 


4.9219 


-1.8947 








6 


4.9877 


4.8986 


-1.7862 








7 


4.9709 


4.8849 


-1.7307 







0.7 


2 


5.6480 


5.5989 


-0.8694 








3 


5.3942 


5.3705 


-0.4394 








4 


5.2684 


5.2503 


-0.3429 






0.5 


1 


5.5967 


5.5510 


-0.8173 








2 


5.2475 


5.2343 


-0.2519 








3 


5.1253 


5.1160 


-0.1821 




0.5 


0.3 


1 


5.3151 


5.2095 


-1.9860 








1.5 


5.2048 


5.1289 


-1.4589 








2 


5.1427 


5.0812 


-1.1966 






0.2 


0.5 


5.2170 


5.1576 


-1.1380 








0.75 


5.1449 


5.0954 


-0.9620 








1 


5.0960 


5.0595 


-0.7163 



Other parameter values 5*0 = 100, r = 0.015, a — 



— 1.6 and T — 1 are used here. 



